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O ' Abstract. 



The results on partial wave analysis of 371 system in reaction Ti^A tvti^ti^A at the momentum 
36.6 GeV/c on the beryllium target are presented. New method of amplitude analysis is suggested 
J> . — extraction of largest eigenvalue of density matrix. Exotic wave with J^'^ = 1 ^p7t is studied in 

O ' four t' regions. No narrow object around M = 1 .6 GeV/c^ is found. Unusually steep t' dependence 

\ for 7t(1300) object is detected. 

(N 

INTRODUCTION. THE PARTIAL WAVE ANALYSIS 

> 

■ The VErtex Spectrometer (VES) setup is a large aperture magnetic spectrometer includ- 
ing the system of proportional and drift chambers, a multichannel threshold Cerenkov 
counter, beam-line Cerenkov counters, a lead-glass y-detector (LGD) and trigger ho- 
doscope. This permits full identification of multi-particle final states. The setup runs on 
the negative particle beam with the momentum of 36.6 GeV/c. The description of the 

^ , setup can be found in [1]. 

■ In this report we present some results of partial wave analysis of the 371 system in the 



reaction % Be ^ 71+71 71 Be for different t' regions \t'\ = 0.01 — 0.07 — 0. 15 — 0.30 — 



D ! 0.80 GeV^/c^. The discussed results are based upon the statistics of about 8.0- 10^ 

^ ' events. Our previous results were published in [2]. 

The PWA has been performed in the 0.8-2.6 GeV/c^ mass region in 50 MeV bins 
^ . for different t' regions. Modified version of the Illinois PWA program [3] with maxi- 

^ I mum likehood method has been used for the analysis. Amplitudes were written using 

isobar model and relativistic co variant helicity formalism according to [4]. Explicit t'- 
dependence f{t) = te~^' was included for waves with nonzero projection of spin on GJ 
z-axis. Density matrix of full rank was used to describe final state. The set of 42 partial 
waves in the form J^LM^ isobar [3] was used in the analysis. Full wave set and param- 
eterization of isobars including special treatment of 7171 S-wave can be found in [2]. For 
the channels with j'^^ = ^, 1++, 2^+ largest waves are assigned to their own density 
matrix elements and are enabled to freely interfere with each other. 



^ Amelin D.V., Dorofeev V.A., Dzhelyadin R.I., Gouz Yu.R, Kachaev I.A., Karyukhin A.N., Khokh- 
lov Yu.A., Konoplyannikov A.K., Konstantinov V.F., Kopikov S.V., Kostyukhin V.V., Kostyukhina I.V., 
Matveev V.D., Nikolaenko V.I., Ostankov A.P., Polyakov B.F., Ryabchikov D.I., Solodkov A.A., Solo- 
vianov O.V., Zaitsev A.M. 
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Wave 2 P0+ p7i at \t'\ < 0.03 GeV^/c^ in a) full density matrix, b) largest eigenvalue; 
pjt at \t'\ < 0.03 GeV^/c^ in c) full density matrix, d) largest eigenvalue. 



EXTRACTION OF LARGEST EIGENVALUE. 

Results of PWA are represented in general by density matrix. For physical analysis 
amplitudes are much more convenient. We present here a new type of amplitude analysis 
— extraction of largest eigenvalue of density matrix. Density matrix can be represented 
by its eigenvalues and eigenvectors: 

4^ w u f ekis k-th eigenvalue, ei> e2> ...> ed 

p= ) ek*Vk*V^ where < t/ ■ i w ■ . 

I ^'^^ eigenvector 

Single out leading term: 

d 

P = pL + Ps, pL = ei*Vi*V{^, ps^Y^ek*Vk*V+ 

k=2 

Here pz, is coherent part of density matrix and p5 is the rest (incoherent part). To be of 
physical meaning, this decomposition must be stable with respect to variations of density 
matrix elements. This is so if eigenvalues are well separated in comparison with errors 
in p matrix: \ei — ^2! ^ ^(Pij)- This is the case for ti+ti^ti" where ei ~ 1, ^2 ~ O-l- 

Extraction of largest eigenvalue has the following advantages. By construction ma- 
trix pi has rank one, so phases are well defined. It quantitatively uses information about 
coherence factors, which is often ignored. Practical experience shows that resonance 
structures tend to concentrate in pi and leakage (see below) is suppressed in pi- Never- 
theless, p5 can contain different non-leading exchanges, albeit it often contains garbage. 
We can also note that if the wave is small in p^, its phase with respect to largest waves 
can not be measured. 

As a restriction this method requires a lot of data for good fit with small errors. It is not 
applicable if eigenvalues are not separated. In this case sometimes a group of clustered 
eigenvalues can be extracted. It is also not applicable if all eigenvalues except one are 
not statistically significant, as it is the case for unnatural sector for Ti'^n^n^ system. 

An example of separation of largest eigenvalue is present in figure 1. On sub-figures 
a) and b) one can see a huge difference between wave 2~P0+p7t in the full density 
matrix and in the largest eigenvalue. It was already noted [5, 2] that at low t' region this 
wave is large and highly incoherent with others at M^^i fti 1.2—1.4 GeV/c^ and have 
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FIGURE 2. 

real data, larj 



Leakage study. Wave 1 +p7t at 0.03 < \t'\ < 1.0 GeV^/c^ in: a) real data, fuU p 
;est eigenvalue; c) leakage, full p matrix; d) leakage, largest eigenvalue. 
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only a relatively small shoulder atM^ji ~ 1.6 GeV/c , which corresponds to well known 
712(1670). Our analysis confirms this as shown in figure 1 a). Contrary to this in the 
largest eigenvalue this wave is dominated by physical 712(1 670). The bump at low M^n 
is at least ten times suppressed, which is consistent with its low coherence with other 
waves. The physical nature of this phenomena is still unknown. 

In figures 1 c), d) one can see the difference between full density matrix and largest 
eigenvalue for the wave 4+Gl+p7i at \t'\ < 0.03 GeV^/c^. At both figures we can see 
at M^Ti ~ 2.0 GeV/c^ a signal for 04 (2050) while in fig.l c) we can also see a bump at 
M3J1 1.3 GeV/c^ which is absent in fig.l d). This bump is a leakage from a2(1320) 
with intensity about 0.2% of total number of events or about 7.5% of 2+ wave in this 
low \t'\ region. This leakage is at least ten times suppressed in coherent part of density 
matrix. 



LEAKAGE STUDY. 

We have used the following method to study a possible leakage effects due to finite 

setup resolution and limited knowledge of setup acceptance. At first, we fit real data 
with small but representative wave set — 12 largest waves. The result of this step is a 
reasonably accurate representation of multidimensional distribution of real events. Next 
we generate Monte-Carlo events according to density matrix from this fit, smear these 
events according to modeled setup resolution and fit them as usual using standard wave 
set with all 42 waves. To study dependence of results of variation of modeled acceptance, 
we have used Monte-Carlo events without smearing, but used corrupted MC program 
(with hodoscope trigger logic excluded) for final fits with 42 waves. 

The results of leakage study for exotic wave 1 '"pji are shown in fig. 2. One can see 
that exotic wave can contain 30-50% of leakage, but can not be described by leakage. 
The structure of density matrix in the real data and leakage is quite different — in the 
coherent part of p leakage is 20-50 times suppressed. Most other waves can contain 
5-10% of leakage. 
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FIGURE 3. Wave 1 Pl"'^pjt at different \t'\ regions: a-d) full density matrix; e-h) largest eigenvalue; 
i-1) phase difference (p(l~+p7t) — (p(2~+/2Jt), full density matrix. 



WAVES WITH JP<^ = 1 +. 

Waves with exotic quantum numbers 7^*^ = 1 +p7l were included in our analysis with all 
possible projections = 1"*", 0~, 1~. Waves withM^ = 0~, 1~ are small in comparison 
with = 1 + (except for Ms^ < 1 .2 GeV/c^) and are not considered by us as significant. 
Wave l~Pl+p7C in four different t' regions is shown in fig. 3. In the full density matrix 
(fig. 3 a-d) at low t' this wave consists mainly of a bump around M — 1.0—1.2 GeV/c^ 
and a shoulder at M Ri 1.6 GeV/c^. At higher t' regions the bump diminishes while the 
shoulder remains approximately the same. From figures 3 e-h) one can see that this 
shoulder corresponds mainly to the coherent part of the density matrix, and this part 
remains stable over investigated t' region. The bins over t' are selected so that numbers 
of events in a2(1320) peak are approximately the same for all bins, namely about 
20000 events/50 MeV, so t' distribution for coherent part of exotic wave is roughly the 
same as for a2(1320) while its crossection is only 2-3% of it. One can see a sharp drop 
on the coherent part of exotic wave at M^-^ 1 .8 GeV/c^. We can see analogous effect in 
some other p% waves and it can be connected with worse description of high M^^j^ region. 
Phase difference (p(l-pl+p(770)7l) -(p(2-S0+/2(1270)7l) is shown in fig. 3 i-1). This 
phase difference is not constant, visible drop corresponds to phase raise of 712(1670) 
resonance. Again the shape of phase variation is more or less stable over inspected t' 
region. In general we can not see here narrow exotic object with M 1.6 GeV/c^. 




FIGURE 4. Wave 7^^ O^+eTi: a) \t'\ = 0.01 - 0.07 GeV^/c^ region; b) \t'\ = 0.15 - 0.30 GeVVc^ 
region; c) relative \t'\ dependence of Jt(1300)/Jt(1800) signals and corresponding Msjt bands. 



WAVES WITH = +. 

Results of PWA for the wave 0~S0"'"£7C in different t' regions are presented in figure 4. 
At low t' two peaks are visible which corresponds to 7r(1300) and 7r(1800). At higher 
t' peak for 7r(1300) is clearly suppressed. Relative strength of 7r(1300) and 7r(1800) 
signals together with relative number of events in M^n = 1.0—1.4 GeV/c^ and M^^ = 
1.6 — 2.0 GeV/c^ bands is shown in fig 4 c) in ten t' regions. One can see that t' 
distribution for 7c(1300) is much more steep than for 7r(1800) or even for the total 
number of events in corresponding M^-^ band. A possible explanation of this phenomena 
can be given if the 71 (1300) bump is at least partially non-resonant and produced by 
Deck-type final state scattering. 

CONCLUSIONS. 

Partial wave analysis of h~^'k~'k~ final state at different t' regions was performed on VES 
data. New type of amplitude analysis was suggested — extraction of the coherent part 
of density matrix. Its advantages and limitations were briefly discussed. 
Wave = l-+p7t was studied in different t' regions t' = 0.010-0.070-0.150- 

0. 300.0.800 GeV^c^. Wave shape is broad and more or less the same in all t regions 
studied. Clear phase variation with respect to 712(1670) is visible in all t' regions. No 
narrow object around M = 1.6 GeV/c is found. 

Waves 7^<^ = 0"+ were studied in the same t' regions. Abnormally steep t' distribution 
for 7u(1300) was established. This phenomena can be understood if J^^ = 0~+ wave in 
7t(1300) region is partially consists of Deck-type background. 

This work is supported in part by INTAS-RFBR 97-02-71017, RFBR 00-02-16555, 
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